Infrared Line Emission in the Interacting Region of Arp 244 (the 
Antennae): Colliding Molecular Cloud Complexes ? 
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ABSTRACT 

We report velocity-resolved spectroscopy of infrared hydrogen recombination lines in the in- 
teracting region of the Antennae galaxies (NGC 4038/4039). It generally has been assumed that 
the active star formation found there is due to the interaction of the disks of the two galaxies and 
indeed two molecular cloud complexes, separated in velocity by ^100 km s _1 , have been observed 
in the southern part of this region. Our measurements imply that the two cloud complexes are 
moving away from each other. This result poses interesting questions about the physical mecha- 
nisms leading to starbursts in Arp 244 and other interacting galaxies. 

Subject headings: galaxies: interactions, kinematics and dynamics, starbursts — infrared: galaxies 



1. INTRODUCTION 

The Antennae (NGC 4038/4039, Arp 244, 
VV245) are a nearby (^22 Mpc) infrared-luminous 
pair of interacting galaxies. Using N-body numer- 
ical simulations Toomre & Toomre (1972) suc- 
ceeded in generating the long tidal tails seen in 
optical/H I images of the galaxies (van der Hulst 
1979). The strongest emission in both the radio 
continuum (Hummel & van der Hulst 1986) and 
the CO (1-0) line observed by Sanders & Mirabel 
(1985); Stanford et al. (1990); Aalto et al. (1995) 
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and by S. W. Lee, K. Y. Lo, R. A Gruendel, and Y. 
Gao (2001, in preparation, hereafter L01) occurs 
in the interaction (overlap) region to the west and 
between the two nuclei, especially in its southern 
dense clouds. The Hubble Space Telescope Ha 
image of the Antennae reveals thousands of gi- 
ant H II regions in the interaction region. X-ray 
observations (Fabbiano & Trinchieri 1983) show 
harder and softer components corresponding to 
SNRs/X-ray binaries and thermal emission from 
gas possibly heated by the interaction, respec- 
tively. These observations all point to there being 
copious numbers of active star formation sites in 
the interacting region. 

Mid- infrared (12—17 /im) observations (Mirabel 
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et al. 1998) of the Antennae by the Infrared Space 
Observatory (ISO) confirm that the starburst is 
taking place largely in portions of the interaction 
region which arc highly optically obscured. The 
ratio of [Ne III] and [Ne II] mid-infrared lines indi- 
cates that stars as massive as 60 M are present 
there (Vigroux et al. 1996; Mirabel et al. 1998). 
The peak of [C II] 158 /zm fine-structure line also 
lies in the interaction region, but overall the ra- 
tio of [C II] to CO is a factor of 2.6 lower than 
usual starburst galaxies, suggesting that global 
star formation is not taking place in the Antennae 
(Nikola et al. 1998). 

Using the BIMA Array with a 7.8" x 6.6" beam 
L01 have found that there are two kincmatically 
distinct molecular cloud complexes separated by 
^100 km s _1 in the southern part of the overlap 
region (see Fig. 1), corresponding roughly to super 
giant molecular complexes 4 and 5 found by Wil- 
son et al. (2000). It seems likely that each of these 
is associated with one of the interacting galaxies. 
Because their relative positions along the line of 
sight are unknown, it is unclear if the complexes 
are approaching or moving away from each other. 
If the redshifted cloud complex is in front, it may 
collide with blueshifted cloud and subsequently a 
starburst, perhaps greater than that currently be- 
ing observed, may occur. If the blueshifted clouds 
are in front, then the complexes are moving away 
from each other and the bulk of their interaction, 
if any, already has taken place in this region. 

Comparison of velocity profiles of properly cho- 
sen emission lines, well separated in wavelength 
but from similarly excited species, could clarify 
the situation, as velocity components associated 
with the two complexes should undergo different 
amounts of extinction. UV and optical lines asso- 
ciated with star formation in these clouds are ex- 
pected to be essentially totally obscured by dust. 
With this in mind we have observed the hydrogen 
infrared recombination lines Pa/3 (5-3, 1.282 /zm), 
Br 7 (7-4, 2.166 fzm), and Bra (5-4, 4.052 izm) from 
part of the overlap region. 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

The observations were made during UT 2000 
28 February-1 March at the United Kingdom In- 
frared Telescope (UKIRT) on Mauna Kea, us- 



ing the facility cooled grating spectrometer CGS4 
(Mountain et al. 1990) equipped with a 256x256 
InSb array. All three nights were clear with sta- 
ble, but not exceptional seeing. We selected the 31 
1/mm echelle, 300 mm focal length camera optics, 
and 2-pixel wide (0.8") slit, which give CGS4 a ve- 
locity range of ^2,000 km s _1 at any wavelength 
in the 1-5 /im region at a resolution of ~16 km s _1 
(two pixels). The ~80" long slit was oriented 
north-south for all observations. Sky measure- 
ments were obtained with the telescope nodded 
600" to the west. 

On the first night, we centered the slit on the 
southern CO peak at J2000 RA = 12:01:55.0, Dec 
= -18:53:01.5, detecting weak line emission at Pa/3 
and Br7, but no emission at Bra. On the follow- 
ing two nights we positioned the slit 3" west of 
the CO peak, which corresponds to the location 
of strongest mid-infrared emission as measured by 
ISO. There we detected all three recombination 
lines at four spatially distinct locations along the 
slit (see Fig. 1), which we call Apertures I-IV (in 
the order of north to south). Details of the loca- 
tions are provided in Table 1. 

The data were wavelength-calibrated to an ac- 
curacy of 5 km s _1 using spectra of arc lamps 
and night sky emission lines and flux-calibrated 
using the A5V star HR 4405 (V=4.10), assumed 
to have the visible-infrared colors given in Toku- 
naga (1999) for such stars. From the point-spread 
functions along the slit we estimate that at the 
wavelength of each recombination line 0.60 ± 0.05 
of the stellar flux entered the slit on each night. 
Spectra were summed over ~10 adjacent rows at 
each of the four spatially distinct line-emitting re- 
gions. Data from the final two nights were com- 
bined to yield the spectra shown in Fig. 2. The 
figure also includes CO spectra obtained by L01 
at those locations. The basic results are summa- 
rized in Table 2. 

3. RESULTS 

3.1. Spatial and Velocity Distributions 

Figure 2 shows broad recombination line emis- 
sion, covering 100-200 km s _1 , at each of the four 
locations along the slit. A faint continuum is de- 
tected in some spectra. Ap I and Ap II are located 
near local maxima of CO emission in the northern 
part of the interacting region, where Wilson et al. 
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12 h l m 56 s 54 s 52 s 50 s 

RA (J2000) 

Fig. 1. — The CO(l-O) moment-zero low resolution (7.8"x6.6") contour map (L01) overlayed on the HST 
3-color (Ha, V-band, and B-band) image. The long vertical crosses indicate the locations along the slit 
where IR recombination lines were detected; from top to bottom these are Apertures I-IV. Solid contours 
represent the red component of the CO emission in the interaction region and the CO emissions in each 
nuclei (contour levels are 6 Jy beam -1 km s -1 times 4, 7, 10, 12, 14, 17, 20, and 25). Dashed contours are 
the blue component of the CO emission in the interaction zone (levels are 5 Jy beam -1 km s -1 times 5, 9, 
13, 17, and 21). 



Table 1 

Parameters of Hydrogen Line-emitting Sightlines 



Label 


Declination Rangc a 


Solid Angle 


A v (from CO) b 






arcsec 2 


Red 


Blue 


Ap 1 


-18:52:07.0 - 18.2 


11.2 


18.0 




Ap II 


-18:52:27.4 - 36.2 


8.8 


15.3 




Ap III 


-18:52:44.8 - 55.2 


10.4 


15.9 


29.0 


Ap IV 


-18:53:01.2 - 09.8 


8.6 


23.9 


22.0 



a At R.A. = 12:01:54.8 (J2000). 

b Visual extinction estimated from CO 1-0 column desnity (see 
text). We use 1530 km s _1 to separate the "Blue" and "Red" 
velocity components. 
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Table 2 
Infrared Line Emission 



Line 




Exposure 




Flux (10"" 


1 W m~ 2 ) a 




Observing 






(sec) 


Ap I 


Ap II 


Ap III 


Ap IV 


dates (2000) 


Pa (1.282pm) 


0.270 


7200 b 


5.1±0.3 


3.9±0.2 


5.3±0.3 


4.6±0.3 


02/29, 03/01 


Br 7 (2.166/im) 


0.115 


4800 


1.06±0.06 


1.08±0.04 


1.72±0.06 


1.41±0.03 


02/29, 03/01 


Br a (4.052/xm) 


0.042 


2600 


4.7±0.8 


4.0±0.5 


10. Oil. 4 


5.8±0.8 


03/01 



"integrated over LSR velocities 1520-1660, 1550-1670, 1320-1495, and 1410-1580 km s \ for Ap I, II, III, and 
IV, respectively 

b Except 6600 for Ap II 



(2000) and L01 have found largely redshifted CO 
emission. Ap III and Ap IV fall into the southern 
part of the interacting region where both groups 
found overlapping red and blue components sepa- 
rated by ~100 km s _1 . 

Fig. 2. — Spectra at Apertures I- IV. The unit of 

flux density for the infrared lines is 10~ 15 Watts 

— ? — i 
m /im . 

In the two northern apertures the velocities 
and profiles of the infrared lines generally match 
those of the CO emission. At Ap I both the 
atomic and molecular line profiles show peaks near 
Vlsr = 1600 km s , close to the velocities of the 
nuclei (the NASA/IPAC Extragalactic Database 
gives radial velocities near 1640 km s _1 for both 
NGC 4038 and NGC 4039), and a blue asymme- 
try in the wings, and at Ap II both show a more 
symmetric line profile peaked near the same ve- 
locity. In the south, however, there are distinct 
differences between the millimeter and infrared 
velocity profiles. Most of the infrared line emis- 
sion at Ap III occurs near Vlsr = 1400 km s _1 , 
~100 km s _1 to the blue of the velocity of CO 
peak emission and in fact further blueshifted than 
almost all of the CO. No infrared line emission 
is present near 1600 km s _1 , the peak of the 
CO redshifted component. At Ap IV the ve- 
locities of peak emission arc ^1500 km s _1 for 
both the infrared and millimeter lines, but little or 
no redshifted emission, which is seen in CO near 
Vlsr — 1600 km s -1 , is detected in any of the 
infrared lines. 

At none of the locations are there more than 
marginally significant differences between the pro- 
files of the three infrared lines. However, the 



signal-to-noise ratios are not high, especially for 
Bra, so modest differences in profiles cannot be 
ruled out. In contrast, the differences between the 
CO and hydrogen profiles are large, especially at 
the southern positions. 

3.2. Extinction 

The CO measurements of L01 and the infrared 
line measurements reported here each provide a 
means for estimating extinction. Using the mean 
CO to H 2 conversion factor determined for our 
Galaxy (e.g., Rand, Lord, & Higdon 1999; Bryant 
& Scoville 1999), 

X ee ^2 = 2.5 x 10 20 H 2 cm- 2 (K km s" 1 )" 1 

: co 

= 6.2 x 10 n H 2 cm- 2 (Jy km s" 1 )" 1 , 

and the gas to dust ratio of Bohlin, Savage, & 
Drake (1978), 

N(H 2 ) = 0.94 x 10 21 A y cm" 2 mag, 

we determine that the visual extinctions through 
the molecular gas at Ap I and Ap II are 18 mag and 
15 mag, respectively. Those at Ap III and Ap IV 
are much larger, Ay <~ 45 mag and Ay ~ 46 mag, 
respectively. All of these estimates incorporate 
both interfcromctric and single dish observations. 
It is known that in some starburst galaxies the 
conversion factor is up to three times smaller 
than in our Galaxy (Smith ct al. 1991; Downes 
& Solomon 1998), and therefore it is possible that 
the extinctions are proportionately less in the An- 
tennae overlap region. 

To derive the extinctions to the infrared line- 
emitting regions, we assumed typical H II region 
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Table 3 
Derived Extinctions 



Line Ratio 




A v 


(mag) 






Ap I 


Ap II 


Ap III 


Ap IV 


Pa/3/Br 7 


1.4±0.5 


3.5±0.5 


4.5±0.5 


4.1±0.5 


Pa/3/Bra 


3.1±0.8 


3.6±0.6 


6.5±0.7 


4.6±0.7 


Br7/Bra 


6.7±2.5 


4. Oil. 6 


10.7±1.9 


5.8±1.6 



conditions (T e = 10 4 K, n e = 10 4 cm 3 ), for 
which the intrinsic ratios are 0.162, 0.0275, and 
0.0777 for Pa/3/H/3, Br 7 /H/3, and Bra/H/3 respec- 
tively (Hummer & Storey 1987). We employed the 
formula obtained by Rieke & Lebofsky (1985) for 
0.9-3.4 fim, along with an interpolation of their L 
and M band extinctions, obtaining factors of 0.270 
at Pa/3, 0.115 at Br7, and 0.042 at Bra relative to 
visual extinction. We used the total integrated 
line fluxes as described in Table 2; the overall gen- 
eral agreement of the three infrared line profiles at 
each aperture suggests that at each location the 
detected lines are attenuated by the same clouds. 

The results are tabulated in Table 3. Signif- 
icantly different values are obtained for different 
pairs of lines at two of the four locations. At 
all four locations the smallest value corresponds 
to the shortest wavelength pair (Pa/3/Br7) and 
the largest value corresponds to the longest wave- 
length pair (Br7/Bra). This indicates that, rather 
than a single source of line emission and an obscur- 
ing screen, there are distributions of line emission 
and attenuating clouds along the lines of sight at 
each location (Puxley 1991). 

4. LOCATIONS OF THE CLOUDS 

The extinctions derived from the hydrogen re- 
combination lines (Table 3) are low compared to 
those derived from the CO measurements (Ta- 
ble 1). This is the case even if the extinctions 
through the molecular clouds are a factor of three 
less (as discussed in section 3.2) except for the 
largest value derived for Ap I. The results indi- 
cate that most of the the observed recombination 
line flux in each of the four apertures is emitted 
near the front surface of the associated cloud com- 
plex and, moreover, for Ap III and Ap IV that 
complex must be the only one that is absorb- 
ing and hence the one that is closer to us. Al- 



though the velocities and profiles of the infrared 
lines at Ap III and Ap IV do not closely match 
those of cither molecular cloud, it is much more 
likely that these blueshifted lines are associated 
with the blueshifted molecular cloud than the red- 
shifted one. Thus we conclude that blueshifted 
molecular cloud is in the foreground. 

A second approach makes use of the lack of red- 
shifted line emission in the region where the clouds 
overlap, although redshiftcd emission is present in 
the northern region where the blueshifted cloud is 
absent. The non-detection at Ap III implies cither 
that no redshiftcd recombination line emission oc- 
curs there or that high extinction makes the in- 
frared lines unobservable. For the latter case, if 
we assume the same ratios of infrared line emis- 
sion to CO line emission in Ap III as in Ap I or 
Ap II we obtain rough lower limits for Ay to the 
redshifted cloud. These may be compared to the 
values of Ay obtained from the CO line strength 
(Table 1) to infer the geometry. We obtain lower 
limits to Ay at redshifted velocities of ^8, ~20, 
and ~27 mag for Pa/3, Br7, and Bra, respectively. 
The most restrictive of these, the Br7 and Bra 
limits, rule out extinction by the redshifted cloud 
alone, even if that emission originates entirely be- 
hind it. We conclude that both clouds would be 
required to produce the extinction and, therefore, 
at Ap III we again infer that the blueshifted cloud 
complex is in front of the redshiftcd complex. The 
above argument could be incorrect if the recom- 
bination line emission from the red-component in 
Ap III is intrinsically weak. At Ap IV it is more 
difficult than at Ap III to clearly separate the two 
cloud complexes based on the CO velocity profile, 
so this test is more problematic than for Ap III, 
although it tends to suggest the same conclusion. 

The significant blueshift of the infrared lines at 
Ap III relative to the blueshifted molecular com- 
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ponent is apparent even at the highest resolution 
(~2") CO measurements of L01. This nearly com- 
plete lack of overlap in velocities is somewhat sur- 
prising. We propose two possible explanations. 1. 
The ionized gas lies close to the front surface of 
the blueshifted molecular cloud complex and is 
largely accelerated outward and toward the ob- 
server. 2. The stars with which the ionized gas 
is associated have a different velocity distribution 
than the molecular cloud. 

Several possibilities exist to test the correctness 
of our conclusions about the relative positions of 
the molecular cloud complexes along the line of 
sight. Measurements of the CO lines at higher 
angular resolution would better match the spatial 
resolution of the infrared measurements. Alterna- 
tively, spectral mapping of the infrared lines would 
allow a better comparison with the present CO ob- 
servations. Searches for redshifted radio recombi- 
nation line emission could determine if any ionized 
gas, undetected at infrared wavelengths, is associ- 
ated with the redshifted molecular gas. 

5. SUMMARY AND CONCLUDING DIS- 
CUSSION 

We have obtained velocity-resolved spectra of 
three infrared hydrogen recombination lines, cov- 
ering a wide range of wavelengths, in the inter- 
action region of the Antennae galaxies, where mil- 
limeter CO observations have revealed two distinct 
clouds of gas with widely different radial velocities. 
Using an empirical extinction formula, standard 
recombination line ratios, and the standard gas- 
to-dust conversion factor, we infer that the two 
molecular clouds currently are moving away from 
each other. This is contrary to the conventional 
expectation of cloud collisions leading to starburst 
activities (e.g. Wilson et al. (2000)). 

The age of the starburst in the interacting re- 
gion has been estimated to be 3-4 Myrs (Hummel 
& van dcr Hulst 1986; Ncff & Ulvestad 2000). If 
this starburst were due to cloud collisions, then at 
present the clouds can only be separated (along 
the line of sight) by 100 km s -1 times a few Myrs, 
or ^300 pc. The cloud complexes arc about one 
kpc in extent in the plane of the sky. If they are 
roughly spherical they would still be physically 
overlapped in the radial direction. Yet we ob- 
serve two well-defined components moving apart 



at 100 km s . How this physical arrangement of 
the molecular clouds might have occurred is not 
clear. There is no easy way to determine the ra- 
dial separation, if any, of the two molecular cloud 
complexes, although in principle, one could model 
the interaction of the gas in Arp 244 in sufficient 
detail to see if this physical situation is plausible. 
The morphology of the molecular cloud complexes 
(Fig. 1) does not suggest physical proximity. 

If the starbursts currently observed are not 
due to direct cloud collisions, the physical mech- 
anisms initiating them require further investiga- 
tion. Current observational evidence of molecular 
gas in starburst regions points to two character- 
istic parameters: high column density (Downes 
& Solomon 1998; Bryant & Scoville 1999) and 
high density (Solomon, Downes, & Radford 1992; 
Downes & Solomon 1998). If high column density 
and high density (both of which are naturally ex- 
pected in the process of gravitational collapse) are 
causal factors for starbursts, then the issue reduces 
to identifying the mechanisms during galaxy in- 
teraction that lead to the concentration of a large 
amount of molecular gas in regions such as the 
interacting region in Arp 244. 

We thank the staff of UKIRT for technical sup- 
port. UKIRT is operated by the Joint Astronomy 
Centre on behalf of the U.K. Particle Physics and 
Astronomy Research Council. This work is sup- 
ported by the National Science Council and the 
Academia Sinica in Taipei. Wc thank the staff 
of ASIAA for helpful discussions and suggestions. 
This work also constitutes part of the Master's de- 
gree thesis of M. C. Liang at the National Tsing 
Hua University. 
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